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Page  116. 

THE  TROPOSPHERIC  REFRACTION  OF  RADIO  WAVES. 

G.  I.  Makarov,  N.  P.  Tikhomirov. 

Introduction. 

This  work  is  the  continuation  of  our  investigations  [1,  2], 
which  concerned  direct  and  reverse/inverse  problems  of  diffraction  on 
the  impedance  sphere  of  a  large  radius.  It  is  dedicated  to  the 
use/application  of  developed  ones  in  [1,  2]  methods  for  the  solution 
of  these  problems  taking  into  account  the  heterogeneous  in  the 
heighr/alrituae  atmosphere.  For  the  investigation  the  exponential 
profile/airfoil,  which  is  a  most  precise  analytical  approximation  of 
the  regular  component  of  the  atmospheric  heterogeneity  of  the  Earth 
[3,  5],  is  selected.  A  question  about  the  role  of  tropospheric 
refraction  was  examined  by  many  authors,  in  this  case  were  studied 
different  profiles/airfoils:  bilinear  [6],  parabolic  [7],  exponential 
[8]  and  the  profiles/airfoils  of  more  general  view,  of  which  was 
required  only  sufficient  "smoothness"  [9-11]. 


The  basic  content  of  the  mentioned  works,  with  the  exception  of 
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[7],  is  confirmation  of  the  validity  of  the  concept  of  equivalent 
radius  [ 9 1 .  The  conclusion  of  the  author  [7]  the  fact  that  for  the 
parabolic  profile/airfoil  the  concept  of  an  equivalent  radius  is 
’■•"org,  erroneous,  as  this  was  noted  in  [14]. 

Page  117. 

The  most  general/most  common/most  total  analytical  results,  which 
relate  to  the  consideration  of  tropospheric  refraction,  are  contained 
in  [11].  The  author  examined  arbitrary  smooth  profile/airfoil,  but 
were  obtained  for  the  eigenvalues  hard-to-visual ize  expressions,  the 
region  of  applicability  of  which  was  limited  by  high  frequencies.  A 
question  about  the  role  of  tropospheric  refraction  at  the  low 
frequencies  up  to  now  remained  uninvestigated.  Furthermore,  both  at 
low  and  high  frequencies  it  is  important  to  have  simple  and 
sufficiently  exact  expressions  for  the  eigenvalues.  This  would  make 
it  possible  to  come  to  light/detect/expose  the  possibility  of  the 
analytical  solution  of  inverse  problem.  The  study  of  these  problems 
composes  the  content  of  this  work. 

SI.  Formulation  of  the  problem. 

We  will  be  interested  in  the  distant  field  of  vertical  electric 
dipole  on  the  impedance  sphere  of  a  large  radius.  Sphere  we  assume  by 
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the  surrounded  inhomogeneous  medium  with  tne  dielectric  constant  e, 
which  depends  only  on  the  radius:  e=e(r).  Subsequently  we  will 
consider  that  the  dependence  e(r)  corresponds  to  the  exponential 
prof ile/airf oil  of  atmospheric  heterogeneity.  Before  extracting  the 
formal  solution  of  problem,  which  it  is  possible  to  construct,  for 
example,  by  the  method  of  normal  waves,  let  us  agree  about  the 
designations:  spherical  coordinates  (r,  br  <p)  is  derived  so  that  the 
surface  of  sphere  coincides  with  the  coordinate  surface. of  r=a,  angle 
-9  is  counted  off  from  the  polar  axis,  carried  out  through  the 
radiating  dipole  (r=b,  0=0).  Emitter  is  characterized  by  the  current  ’ 
of  amplitude  I,  which  changes  in  the  time  according  to  harmonic  law 
e-"’u,  and  by  effective  height  ''i  Let  us  extract  now  expression  for  the 
vertical  component  of  electric  field \Er  in  the  form  of  the  series/row 
of  the  normal  waves 


r  l*vJ'hK  |  y» 

\(\tl)(2VM) 

'  4  »' «*(*>*«.<'’>  ~ J 

-  - 

*  -u 

a 

Or 

O'i 

L~frT  J 

V—  V 

5 

/?v  (bMt  (r)  r,  [com- -oil 

_ J _  S  _ l _ 

'  R;  (“)  M"  V 


0) 


Here  relative  dielectric  constant  of  the  atmosphere? 

e 0=8 . 85 *10 ‘ 1 2  F/m,  Mo=4m*10*'’  K/m;  k0=cji/e0Me  -  wave  number  in  the 
vacuum,  p\  |cos'- -  0)|  _  Legendre's  function:  ^(r)—  eigenfunctions;  — 
the  eigenvalues  of  radial  operator  which  corresponds  to  the 
boundary-value  problem 


1)1^,  (0  =  0. 


(2) 
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1)  -jr  =  —  ik.fi’R" J'O,  -)  R.  (r)  £ l.2 1«,  co)  when  lm  >),  where 

ur  k  \  ' 


.•  ,  i  ,o.n  I  .  given  surface  impedance  of  sphere;  *»>  (")■ 


relative  dielectric  cor....  ant  on  its  surface. 


Differeitial  expression  l-  takes  the  form 


L'  r2  {!*  +  ^r)  ~  V  *« 


Con.'  trine  it  carried  out  of  the  condition 


V,i'N  1  1  Vj|(K - G">  1,  [til  Vv  (rr  —  0)>1, 

let  us  lea^  ^1)  to  th-.  standard  form 


_  i-  <  ' 

r  *w/v  1‘hff  i  f  2  i  /"  0  o 

E'~  -  --r-1  •  v  >  7iF«  *  **■/.<*>/,  (r)  (4) 


Are  here  introduced  the  following  designations: 


excitation  coefficient  s— of  the  normal  wave 


vji?(v,+  1)(2',  +  1) 


a(ka)™±  \r-a 


the  high-altitude  factors 


\rt~Vi)  ,i 2  /,f-  (h) 

f  <b)~  *  m  ’  ’  •  '  - 


/  / ^  .  V *m  (")  ^n(r) 

7,1  VC<o  '  r2  (<I)- 


The  important  characteristics  of  normal  waves  is  their  phase 


speed  along  the  surface  of  sphere  and  attenuation  factor  5 ' ■ 
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For  the  s-th  mode  we  have 

< «»  /-''I 

1 "  =  Hi-  <  •  (0) 

'ii _ I  in 

7  ~~~  («) 

It  is  known  that  the  solution  of  transcendental  equation 
relative  to  the  marks  of  addition  v«.  Composes  basic  difficulty  during 
further  investigation  of  the  series/row  of  normal  waves  (4)  In  the 
present  work  it  is  analogous  how  this  was  done  in  [1],  we  will 
attempt  to  obtain  simple  analytical  dependences  v,  on  impedance 
frequencies  w  and  the  parameters  of  the  heterogeneity  of  the  <V 
atmosphere.  In  contrast  to  [1]  we  will  be  interested  only  in  the 
first  approximation  of  the  Galerkin  method. 

Page  119. 

This  method  is  applicable  upon  consideration  of  inhomogeneous  medium 
since  under  the  impedance  conditions  on  the  surface  of  sphere  the 
spectrum  of  operator  lr  is  discrete/digital-.  It  is  assumed  that  the 
spectrum  of  simple.  For  the  realization  of  the  Galerkin  method  more 
conveniently  to  switch  over  from  to  integral  operator  which  proves 
to  be  limited.  Before  carrying  out  this  transition,  let  us  rewrite 
(2),  conditions  1)  2)  in  the  new  dimensionless  variable 

x  ~.k(r  —  a). 

where  k=~-kn and  (0)  =  sm (a*) !r_„—  dielectric  constant  on  the 

surface  of  the  sphere: 
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I  <P  ’m  '  •'  ) 
_Mj<  -f-J) 


\  s„(-v) 


<l‘  1 _ 

,/'3  V  m  U) 


-|  /?,4(.v)  =  0. 


(7) 


.  rfff, 
1)  .* 
d  i 


=  -/o/^(°)>  2)t 


when  fm  k  >  0  (*)  £  ij  (o, 


Here\.~  modified  impedance,  which  depends,  as  is  evident,  from 
the  value  of  the  gradient  of  f unct ion tm  (x)  y on  the  surface  of  sphere, 
equal  to 


(8). 


§2.  Transition  from  operator  Lr  to  the  integral  equation. 


Analogous  to  [1]  we  consider  the  relationship/ratio 


satisfied  for  all  x.  For  the  case  of  homogeneous  medium  the  use  of 
tais  relationship/ratio  in  (7)  corresponds  to  transition  from  the 
cylindrical  functions  to  their  asymptotic  representations  through  the 
Airy's  functions.  The  possibility  of  this  transition  for  the 
examined/considered  case  ka>>l  is  proved  in  the  work  of  V.  A.  Fok 
[12],  Introducing  certain  initial  approximation/approach  v^’  to 
unknown  eigenvalue  v*and  taking  into  account  (9),  let  us  write 
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equation  (7)  in  the  form 

/-A  (x)  ~h  l-Jl  ^  (■<)  ~  0, 


where 


/  ,  ,fl  4. 1  .!  L"< ((1). 0  /.  2.*_\ 


1  (0)  ^  ~nt(x)  , - — — 

,a  V£m(.r) 


(*«)2 


('-£)■ 


(10) 

(11) 


(12) 


Page  120. 


To  differential  expression  (11)  is  related  the  linear  part  of 
function  ,  moreover  before  rfW'r  *s  set  SY11^01  5 ^  |  ^  •  Symbol  6  is 
introduced  with  that  target  in  order  not  to  write  two  different 

cm  (I  ) 

equations  for  the  cases,  when  to  Lx  one  should  carry  x  and  when 

this  to  make  inexpediently.  For  the  reduction  of  recording  (11)  (12) 
it  is  convenient  to  introduce  the  designations 

o  _ -M--t  t- 1)  ,/  .  -((  •'!  ‘) 

■ '  ''r  (kaY  ’  ‘ 5  • '  Yn2 


>>£)/ 

—  I  XT  +  *  • 


-P. 


ka  1  ' 

— 7,  =  a’ (l  —  p{).  f  s  i  (s’) 


1 7 
0„ 


and  the  new  variable 

_  *'2 

v=  (aj)  x  +  ’is-  (13) 

After  this  we  turn  operator  Lir  by  which  we  understand  the  operator, 
formed  by  differential  expression  (11)  and  boundary  conditions  1)*, 
2)*,  and  we  obtain  the  integral  equation 

At  (y)  +  ^  [  AR,s  (0  +  CR,t  (0  j  +  < 0  =■  0-  (14) 
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Through  A,  C  and  B  the  integral  operators 


AR,t(y)~  j/v(v.  f)R,  (t)<ift 

I* 

v 

CR.t  (y)  s  -  ~  (cl)1'2  j‘  A'(y.  /)  (/  -  T,)  /?.,  (0^. 

*  c 

BR,s(y)~  |  K  (y,  /)  «l>  (/)■/?.,  (/)<//. 


are  designated.  Furthermore,  the  nucleus 

AT  (y,  t)--~  1  'My),M»(/).  7,  <  t-  v  (is) 

ju,  (y)-«,(/),  /, 

«i (O’)  =  -^7  1°  «,(v)-//2(y)i,  (19) 

a  ^  'MtO  +  •!•«* t  (20) 

«,  rrJ  -1-  </*"i  (y,) 

ua  and  Uj  -  two  linearly  independent  solutions  of  the  equation 

%+y»=o.  (2D 

moreover  Ki (y) G It,.  I.  and  u2  (y )=~u1  (y ) . 


Page  121. 


Function  F(t),  entering  in  (17),  is  determined  by  the  equality 


_ Li/=  (n  'l2- 

i  y  m  y*-/  (4*i 


Instead  of  the  duct/contour  of  integration  l.<  (Fig.  1)  let  us 
introduce  the  equivalent  to  it  F,  (Fig.  2),  on  which  function  ua(t) 
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with  t-*°=  exponentially  vanishes,  providing  the  .limitedness  of 
operators  A,  C  and  B  and,  consequently,  also  the  applicability  of  the 
Galerk?n  method  for  finding  the  eigenvalues  aX  of  equation  (14). 

§3.  General/common/total  expressions  for  the  eigenvalues. 

Without  defining  concretely  thus  far  the  form  of  the  function 
will  obtain  general/common/total  expressions  for  eigenvalues 
v.  in  the  first  approximation.  We  search  for  v<in  the  form 

v,  =  fci(l  +  £\  (23) 

and  in  the  same  form  we  assign  the  initial  approximations/approaches 

v(  =  ku  ^  1  + 

Subordinating  to  condition  a’>l.-we  will  have 

■X  <-r>> 

Setting  for  itself  as  a  goal  to  obtain  sufficiently  precise  and 
simple  analytical  dependences  >,  on  the  parameters  of  the  problem,  let 

us  examine  only  two  limiting  cases  -  (I)-,  when  and  (II)-,  when 

\(],  >1 
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Without  repeating  those  given  in  [1]  reasonings  relative  to  the 
selection  of  initial  approximations/approaches  t{  (j=0)  for  the  case 
(I),  j=®  for  the  case  (II),  let  us  recall  their  definition: 


"I  ( —  0  =  °-  (26) 

(5  =  0,  1,  2,  ...) 

«,(  —  /;)  =  o.  (27) 

As  the  coordinate  element/cell  z0  we  select  function  ul(y),  which 
satisfies  condition  (26)  in  the  case  (I)  or  (27)  -  in  (II).  Then  in 
the  first  approximation,  for  the  eigenvalues  we  obtain 

t  —t>  =  >1-  "ill  (jEh  IA  .nov 

'  '  (^,.5i) +  (<:«,.  at)  ■  (28) 

Let  us  turn  to  formulas  (18) -(20),  that  determines  nucleus  K(y,  t). 
It  is  not  difficult  to  note  that  a  is  the  high  parameter  both  in  th« 
first  and  in  the  second  case.  Actually/really,  when  we  have 
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Then  formula  (28)  can  be  written  thus: 


'■  1,- . 

(32) 

Here 

r 

(33) 

V 

j  (t  —  I<)  (0  dt  — ( *1) 

s 

(•'.-}) 

V-  *tj  <I>  (0  xi\  (/)  (It 
r* 

(35) 

Page 

123. 

From  the  designations 
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introduced  earlier  it  is  evident  that  under  condition  a’\M.  which  is 
considered  carried  out,  s  'Taking  into  account  this,  expression  for 


let  us  wr?te  in  the  form 


(36) 


where  ,  ,  - (0)  ka  equivalent  radius,  which  is  eoual  to  a  radius 

of  sphere  a  for  the  homogeneous  medium,  when  sm(0)=0.  in  this  case 

)  -Relatively  we  also  assume  that  After  returning  to 

formulas  ( 32 ) - ( 34 ) f  we  see  that 

l+£-=<+T(vf  i-  <37» 

Second  term  in  (37)  is  negligibly  small  in  comparison  with  one  for 
several  first  numbers  s  by  force  V>  I,  therefore  subsequently  we  will 
assume  that 


'  +  %  = 


(38) 


Formula  (32)  taking  into  account  (29)  (30)  and  (38)  is  converted  to 
the  form  for  (I)  case 


(39) 

(40) 


During  the  conclusion/output  of  these  relationships/ratios  it  is 
taken  into  consideration  also  that  factors  x0  and  ya>*  which  appear 
before  first  terms  in  (39)  and  (40),  are  in  effect  equal  to  one  [1], 
It  is  obvious,  f irst  term  in  (39)  (40)  corresponds  to  the  equivalent 

n ' ' 

"uniform"  problem,  in  which  instead  of  radius  a  and  impedance  \ 
figure  equivalent  radius  «3  and  modified  impedance  V 
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Component/term/addend  ~  in  (39)  (40),  generally  speaking,  appears  as 
a  result  of  a  difference  in  prof i  1  e/airfoil  ‘■‘-(•'O  from  the  linear; 
however,  it  is  not  equal  to  zero  even  for  the  linear  profile/airfoil. 
The  estimation  of  this  component/term/addend  with  linear  function 
*-'»i(*)  is  conducted  in  the  following  paragraph.  At  conclusion  of  this 
paragraph  let  us  note  that  if  we  are  interested  in  the  wider  range  of 
chance  in  parameter  q\  then  instead  of  first  term  in  (39)  (40)  it  is 
possible  to  substitute  the  appropriate  expressions  from  [2],  obtained 
in  the  second  approximation/approach  of  the  method  of  moments/torques 
for  the  homogeneous  medium.  For  this  should  be  used  formulas  (8)  and 
(9)  from  (2),  having  preliminarily  replaced  in  them  q  by  q». 

Page  124. 

§4.  Linear  profile/airfoil. 

For  the  linear  profile/airfoil  the  proof  of  the  concept  of  an 
equivalent  radius  usually  is  conducted  at  the  level  of  differential 
equation  (7).'  In  this  case  the  assumption  about  the  possibility  to 

^  /  — . —  d-  '  ' 

disregard/neglect  the  so-called  derivative  of  Schwarz  V  £m (■*> dxi  • 

Substantially  is  utilized* It  is  considered  also  that  upon  transfer 
from  the  problem  with  the  unhomogeneous  atmosphere  to  the  the 
equivalent  to  "uniform"  the  impedance  of  sphere  does  not  change.  In 
actuality  in  the  equivalent  "uniform"  problem  must  figure  the 
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modified  impedance,  which  can  considerably  differ  from »V  by  very  low 
frequencies.  It  is  of  interest  to  rate/estimate  the  effect  of  the 
named  factors  on  eigenvalues  /,•  Since  for  both  cases  (I)  and  (II)  in 
question  in  the  distant  zone  the  field  is  determined  by  one  normal 
wave,  to  ’..nich  eigenvalue  t„  corresponds,  we  will  obtain  analytical 
expressions  only  for  this  eigenvalue.  For  the  same  reason  and  for  the 
exponential  profile/airfoil  we  will  be  interested  only  in  zero 
eigenvalue.  Considering  that 

(x )  =  sm  (0)  4-  (0)  x,  (■!’.' 

entire  function  e,n(x)  is  related  to  operator  (11),  for  which  we 
assume/set  5=1.  In  this  case  we  have 


(42) 


For  future  reference  it  is  convenient  to  introduce  parameter  rj0 
according  to  the  formula 


n* »  — 


«!>> 
(*  ) 


(43) 


We  will  consider  it  its  small  in  comparison  with  one:  M<l.  Under  the 
conditions,  characteristic  for  the  earth's  atmosphere,  this 
parameter,  for  example,  has  a  value  -0.0063  at  the  frequency  f  =  5  kHz 
and  with  an  increase  in  the  frequency  decreases,  vanishing  with  f -»«*>. 
Let  us  examine  now  coefficient  of  b„,  which  determines  second  term  in 
(39)  (40)  with  s=0 

b„=  j  <D  (t)uU  )dt  (44) 

T. 
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Page  125. 


Here 


<I»(n=r 


4*  <;  |  |  -  r„  (/  -  T,„p  • 


(45) 


Let  us  recall  that  T<. ~  —  t’a{j  —  0,  '«-).  to  express  integral  (44)  through 
the  known  functions  is  impossible,  but  it  is  possible  to  calculate 
approximately  it.  For  this  we  will  use  smallness  tfo  and  fact  that  the 
basic  contribution  to  the  integral  gives  region  small  t.  This  follows 
from  the  character  of  the  behavior  of  function  U! (t) ,  qualitatively 
depicted  in  Fig.  3,  where  X  -  certain  complex  coefficient. 
Decomposing/expanding  F(t)  in  the  series/row  according  to  degrees  ^ 
it  is  not  difficult  to  obtain  correction  to  eigenvalue  t0  due  to 
Schwarz's  derivative.  We  v-ill  be  restricted  to  the  dominant 
relatively  tfo  term 

"  — -XYi’;[l  +0  (^)l-  (46) 


Let  us  turn  now  to  the  expression  for  Taking  into  account  (8)  and 
(43),  let  us  write  it  in  the  form 
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As  is  evident,  </’  in  the  accuracy  corresponds  to  the  concept  of  an 
equivalent  radius,  since  it  is  obtained  from  q  for  the  homogeneous 
medium,  replacement  of  radius  a  on  Component/term/addend  in 

(47)  is  correction  to  the  concept  of  an  equivalent  radius  because  in 
the  presence  of  inhomogeneous  medium  instead  of  impedance  appears 
modified  impedance 
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Fig.  3. 

Page  126. 


It  is  clear  that  this  component/term/addend  can  prove  to  be  essential 
only  in  small  ones  We  will  obtain  the  resultant  expressions  for 
t0,  after  substituting  (46)  (47)  in  (39)  (40)  with  s=0 


(i  )h -  C  ~  £  -h  ~i  -n.  -  -r  x !  i  +  o  (’i.)l. 


(ID  t0-tZ 


-  + 


(19) 

m 


Thus,  even  with  the  linear  profile/airfoil  the  concept  of  an 
equivalent  radius  can  be  considered  justified  only  in  such  a  case, 
when  in  (49)  (50)  the  terms,  which  contain  parameter  are 
sufficiently  small.  Let  us  rate/estimate  their  effect  on  the  phase 
speed  and  attenuation  of  zero  normal  wave,  for  this  let  us  examine 
the  linear  profile/airfoil,  characteristic  for  the  lower  layers  of 
the  earth's  atmosphere 


-V 


DOC  =  84095300 


PAGE  19 


with  a=6.4*10"4,  H=7.4  km. 


Latter/last  component/term/addend  in  (49)  (50),  that  appears  due 
to  Schwarz's  derivative,  obviously,  does  not  affect  fading;  it  it  is 
possible  not  to  consider,  also,  during  the  calculation  of  phase 
speed,  since  it  the  second  order  of  smallness  in  comparison  with  the 
contribution  due  to  the  modified  impedance.  For  the  certainty  we 
assume/set  — 0.  Table  1  gives  values  and  af°yk0  without  the 

account  and  taking  into  account  second  term  in  (49)  for  several 
frequencies . 
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'Table  1. 


Note;  a)  without  the  account;  b)  taking  into  account  second  term 
in  (49)  for  several  frequencies. 

Key:  (1) .  kHz.  (2) .  s. 

§5.  Exponential  profile/airfoil. 

For  the  exponential  profile/airfoil  we  will  examine  the  cases  of 
hicn  and  low  frequencies  separately.  The  reason  for  this  separation 
consists  in  the  fact  that  eigenvalues  for  these  cases  substantially 
differ  from  each  other.  It  proves  to  be,  and  this  it  was  possible  to 
assume,  on  the  basis  of  the  physical  considerations  that  at  the  high 
frequencies  the  eigenvalues  for  the  exponential  profile/airfoil  were 
close  to  the  appropriate  eigenvalues  for  the  linear  profile/airfoil, 
and  on  the  low  frequencies  they  were  located  near  the  eigenvalues  of 
problem  with  the  homogeneous  medium. 


Page  127. 


rvTvwrf\>  v-"vvvv' v*  v*  .-»  v* *.%  ■*  \  /.  % /.  .*.  %  ,*.  *.  /. 
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Therefore  it  is  necessary  to  introduce  different  initial 
approximations/approaches  v  for  the  high  and  low  frequencies  in  orde 
already  in  the  first  approximation,  of  the  Galerkin  method  to  obtain 
sufficiently  precise  formulas  for  unknown  eigenvalues  w.  As  already 
mentioned  in  the  preceding/previous  paragraph,  we  will  be  restricted 
to  the  investigation  of  zero  eigenvalue  v0,  moreover  let  us 
disregard/neglect  Schwarz's  derivative.  Based  on  the  example  of 
linear  profile/airfoil  it  was  shown  that  its  contribution  to  the 
eigenvalue  of  second  order  of  smallness  in  comparison  with  the 
contribution  of  the  modified  impedance. 


High  frequencies.  For  the  high  frequencies  the  linear  part  of 
function  given  by  the  formula 

r 

sm(x)=\-\-*-e~UT ,  (52) 

it  is  carried  to  operator  L,(ll),  assuming/setting  5=1.  Coefficient 
b0  in  this  case  determines  correction  to  eigenvalue  t0  due  to  a 
difference  in  the  exponential  profile/airfoil  from  the  linear  and 


takes  the  form 


f  ! - 1  +  Mf -  To)}  K? (0 dt.  (53) 


DOC  =  84095300 


PAGE  22 


Equivalent  radius  «•  is  determined  by  formula  (42).  We  will  consider 
high  those  frequencies,  with  which  is  fulfilled  the  inequality 

i. 


In  order  to  represent  the  real  values  of  parameter  is  let  us  note  that 
under  the  conditions  of  the  earth's  atmosphere,  when  a-6.4*10''1, 

H-7 . 4  km,  we  have  |t,  o,Of>  at  frequency  f=10  MHz.  With  an  increase  in 
frequency  is  it  decreases,  vanishing  with  f-»®.  Three  times  integrating 
in  parts  first  term  in  (53)  and  utilizing  smallness  of  parameter  is.  it 
is  not  difficult  to  obtain  the  approximation 


—  ~  Y r  -  --2/7  !S  5  (*fo) “H  Vh  (m>)]  +  10  ui(lo)j  X 

X  1 1  o  (is)l-  (54) 


For  conclusion/output-  (5*)  recursion  relations  (ri8)-(ni0)  of  work  [1] 
were  used,  and  it  is  taken  into  consideration,  that  ux (7, ) *u' 2 (70 )=0 
in  view  of  the  selection  of  initial  approximation/approach.  It  is 
obvious,  with  cj-*®  coefficient  bo-*0,  therefore,  with  u-*°>  disappears 
the  difference  between  the  exponential  and  linear  profile/airfoil. 
This  bears  out  the  fact  that  only  the  lower  layers  of  the  atmosphere, 
where  the  prof ile/cirf oil  is  close  to  the  linear,  rlay  main  role  in 
shaping  of  the  distant  field  of  high  frequency. 


Page  128. 

Substituting  (54)  in  (39)  (40)  and  taking  into  account  (47)  (33),  we 
come  to  the  final  formulas  for  eigenvalue  t0  in  the  case  of  the  high 
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frequencies 


,0V! 


(i)  t0  —  -(-  — ;j  Tj,  -f-  - 

^*ft 


3-.r)  (*<•) 


i 

TTT 


1  +  T  2  H 


lSll  +0(^)1,  (55) 


do  4-^^^  +  r^---^~)(/(r)V,ii+o(il,)!.  (56) 


Table  2  gives  formulas  (55)  obtained  with  the  use  of  the  value  of 
phase  speed  and  attenuation  of  zero  normal  wave  when  t\  =  0  for  the 
linear  and  exponential  profiles/airfoils. 


Low  frequencies.  Assuming  that  at  the  low  frequencies  the  role 
of  the  atmospheric  heterogeneity  of  exponential  profile/airfoil  is 
not  so/sn'-h  essential  as  on  the  high  ones,  we  will  search  for 
eigenvalues  in  the  vicinity  of  the  eigenvalues  of  problem  with  the 
homogeneous  medium,  in  which  k=k0.  For  this  of  the  function  p„,{x)  to 
operator  Lj  (11)  one  is  carried,  assuming/setting  5=0,  and  by  k 
everywhere  we  understand  k0.  Then  coefficient  b0,  which  determines 
the  effect  of  tropospheric  refraction,  is  written  as  follows: 

bn  —  (ni  j  e"’u-"\i\{l)dt.  (57) 


Here 


(¥f 


Low  we  will  call  those  frequencies,  with  which  is  fulfilled  the 
inequality 
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For  the  earth's  atmosphere  at  the  frequency  f=5  kHz  we  have  m=8.94. 
With  the  decrease  of  frequency  the  parameter  u  increases.  Three  times 
integrating  in  parts  (57)  and  considering  as  that  carried  out 
relationship/ratio  (58),  easily  we  come  to  the  approximation 


h _  :l1  I  ,  o  Mi "!'<•'>) 

" |~ - r  - - 


-I-  0  ( ~rjj  ■  (W 
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Table  2. 


20  0.021002  0  021741  78,41  78,05 

50  0,016008  0.0 1 0064  106,47  106.84 


Note:  a)  for  the  linear  profile/airfoil;  b)  for  the  exponential 
profile/airfoil . 

Key:  (1) .  MHz .  (2) .  s. 

Page  129. 


It  is  obvious,  with  the  decrease  of  the  frequency,  when  m-*00, 
coefficient  bo-»0,  and  heterogeneity  of  the  atmosphere  in  this  case 
becomes  apparent  only  through  the  modified  impedance.  After 
substitution  (59)  in  (39)  (40)  we  obtain  the  approximation  analytical 
formulas  for  t0/  valid  at  the  low  frequencies 


(I)  ■  +  4r’l  +  «. 


(00) 

(61) 


During  writing  (60)  it  is  taken  into  consideration,  that  in  this  case 
=  since  —  «0.  In  formulas  (60)  (61)  term  o(l/#i)  is  small  in 

comparison  with  one  not  only  because  ^>>1,  but  also  due  to  the 
character  of  the  behavior  of  integrand  in  the  rejected  integral.  This 
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provides  applicability  (60)  (61)  in  the  wider  frequency  region,  when 
n>  1 . 


§6.  Solution  of  inverse  problem. 


In  our  work  [2]  inverse  problem  for  the  impedance  sphere,  which 
is  located  in  the  homogeneous  medium,  was  solved  by  amplitude  method. 
The  use/application  of  this  method  for  solving  the  inverse  problem 
taking  into  account  the  atmospheric  heterogeneity  of  exponential 
profile/airfoil  composes  the  content  of  present  paragraph.  We 
consider  source  and  receiver  as  those  arranged/located  on  the  surface 
of  sphere  and  as  in  [2],  we  will  solve  problem  under  the  assumption 
of  single-modality  of  the  propagation,  when 


A-iOfc,  ,„•«  *a,A  ,v,„ 

n'  4 a  V  nfaiQ  V  MllO*  e  • 


(62) 


Here  x  -  given  distance  from  the  source  to  observation  point 


ka  1 
1 


T-ftTO’*- 


(63) 


Taking  the  ratio  of  amplitudes  |£',(jc,)!  and  l (-v2) !.  measured  at 
frequencies  <»<  at  two  diverse  points  x3  and  x3,  we  will  obtain  system 
of  equations  for  the  determination  of  the  electrical  parameters, 
which  determine  impedance,  and  the  parameters  of  the  heterogeneity  of 
the  atmosphere 

f-'-f  ’  Im  tn  | p„,  »f|  =  •!»,, (/,  A  =  I,  2,  . . . ,  /V),  (64) 
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where 


l-r  (-'l ) 


(63) 


and  pk—  unknown  parameters,  whose  number  is  equal  to  N. 


Page  130. 

With  known  analytical  expressions  /of/?/. .  «,•]  the  solution  of  system  (64) 
will  not  compose  labor/work  and  can  be  obtained,  for  example,  by 
Newton-  Kantorovich  method  [13].  The  necessary  expressions  for  t0  are 
obtained  in  the  preceding/previous  paragraph  in  two  limiting  cases 
(I)  and  (II)  both  or.  the  high  ones  and  at  the  low  frequencies.  If 
necessary  it  is  possible  to  utilize  also  formulas  (8)  (9)  work  [2], 
as  already  mentioned  earlier. 


Here  we  will  pause  it  the  case,  when  system  (64)  admits 
analytical  solution  relative  to  parameters  a  and  H.  let  us  examine 
the  frequency  band,  bounded  above  by  conditional,  and  from  below  by 
the  frequencies,  at  which  still  it  is  possible  to  assume/set  so 

that 


'o- 


— 


i  -t  h- 


'h-  —  V 


(66) 


In  the  lower  part  of  this  range  with  u»l  the  effect  of  tropospheric 
refraction  is  negligibly  small;  therefore  in  this  part  of  the 
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V  -**’ «w 


V-  .v ^ 
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spectrum  it  is  possible  to  determine  the  impedance  of  sphere  *\ 


counting  the  medium  of  uniform.  This  problem  is  solved  in  [2]; 


therefore  we  will  consider  impedance  as  the  function  of  frequency 


known.  However,  for  determining  a  and  H  we  will  use  the  higher 


frequencies  of  the  range  in  question,  when  the  effect  of  tropospheric 


refraction  is  noticeable.  Taking  into  account  (66)  from  (64)  we 


obtain  system  of  two  equations  relatively  a  ana  H 


a  2//*  —  ?  ^  ~  1  »  2)- 


Function  <pK)  is  determined  by  the  equality 


?  <->  (*<) + "»  ('2+ -£))  <68> 

«Jn  |  /  k^a  ,1/3 

and  it  is  considered  known.  Parameter  i(u>) s =  ^ I  — )  •  The  solution 


of  system  (67)  is  located  elementarily,  and  we  immediately  extract  it 


f  ?l  _  Jt_ 

tj _ _j!_  _ n_ _ 

2,n 


2  H*  +  i\ 

7.  - '0  • 

«1?  •" 


For  the  reduction  of  recording  the  designations 


T  (-»/)  =  Tf.  ?("><)  =  ?/• 


are  here  introduced.  In  formula  (70)  the  selection  of  mark  i, 


naturally,  does  not  play  role. 


Page  131. 


/  ~_r  V  *  *■ 


■’v. '»V v‘. 


•  v  *>  V''-^  \v’>  v*wi>  \> 


.  •  tv  ;  ,s(  i"v ".s.  V-  '.-v '. '•  gsi  --v  >.•>  £'*  .•>  ; 
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Page  141. 

REFINEMENT  OF  BOOKER'S  EQUATION  FOR  CYLINDRICAL  LAYERED  INHOMOGENEOUS 
ANISOTROPIC  MEDIA. 

V.  M.  Vyatkin. 

The  known  equation  of  Booker  [1]  is  obtained  for  the  case  of 
plane-layered  anisotropic  medium.  But  the  real  ionosphere  possesses 
final  curvature  and  terrestrial  magnetic  field  is  heterogeneous. 
Therefore  it  is  of  interest  to  refine  Booker’s  equation  taking  into 
account  these  two  effects.  We  selected  the  cylindrical  model  of  the 
ionosphere,  which  makes  it  possible  in  the  first  approximation,  to 
consider  the  curvature  of  the  real  ionosphere. 

Let  us  introduce  cylindrical  coordinate  system  (p,  cp,  z)  and  we 
will  consider  the  ionosphere  as  layered  inhomogeneous  on  p  the 
anisotropic  medium,  which  occupies  the  region  of  space  with  p» b.  We 
consider  suosequently  earth's  magnetic  field  heterogeneous,  which 
depend  only  on  p  (heterogeneity  of  field  in  the  directions,  the 
tangents  to  the  earth's  surface,  as  are  shown  calculations,  can  be 
disregarded/neglected) . 
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The  general/common/totai  equation 

rot  ro:  E  —  mE  =  0,  (1) 

where  k0=co/e0M0  -  the  wave  number  of  vacuum,  is  initial  in  the  study 
of  the  problem  about  the  propagation  of  electromagnetic  waves  in  the 
anisotropic  media,  t.,;—  tensor  of  the  relative  complex  dielectric 
constant  of  medium. 

The  solution  of  equation  (1)  we  will  search  for  in  the  form  of 
series/row  according  to  the  reverse/inverse  degrees  of  k0. 
Substituting  this  series/row  in  equation  (1)  and  equalizing  the 
coefficients,  which  stand  in  the  equal  degrees  of  k„,  we  obtain  the 
chain/network  of  the  engaged  equations,  the  first  of  which  takes  the 
form 


[V'J-jvT-E 


«or 


Ji 


■VmE'°'  =  0. 


(2) 


It  is  not  difficult  to  show  that  in  the  case  of  cylindrical 
layered  inhomogeneous  (on  p)  anisotropic  medium  interesting  us  the 
vector  M  has  components 


vT  —  ;  -7  C:\  C2 1, 


(3) 


where  Cj=cos  7„  and  C2=cos  a,  and  a  has  the  sense  of  the  angle 
between  the  direction  of  wave  and  z  axis;  ^ „  can  be  considered  as  the 
angle,  formed  by  zonal  harmonic  [2]  (by  spiral  wave)  of  order  v  with 
the  direction,  by  tangent  to  the  cylinder  of  radius  k0p0(cos 
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7„  =  j>/k oPo);  P o  -  coordinate  of  the  point  of  the  position  of  the 
source,  which  we  will  assume/set  by  that  arranged/located  lower  than 
the  boundary  of  ionosphere  (p0<b). 

Page  142. 

For  value  94,/3p=q(p)  from  the  condition  for  existence  of  the 
nontrivial  solution  of  equation  (2)  the  equation,  which  is  the  analog 
of  Booker's  equation  [1]  in  the  cylindrical  case 

F {<])  =  aq'  -f  pi/3  -(if2  -|_  ?xf  -|_  e  =  0.  (4  j 

is  obtained. 

Let  us  note  that  equation  (4)  is  obtained  from  "Cartesian" 
biquadraticof  Booker  [1]  by  replacement  of  Cx  by  (p0/p)C , ,  which 
follows  from  comparison  (3)  with  the  expression  for  components  v'F  in 
the  flat/plane  case.  Therefore  expressions  for  the  coefficients  a,  /3, 
7,  5,  e  of  equation  (4)  take  the  same  form,  as  in  flat/plane  case  [1] 
with  replacement  of  Cx  by  (p0/p) C,. 

If  we  in  the  formulas  for  the  coefficients  of  equation  (4) 
assume  y  ^  1  _7i.  where  7?  =  5/a  (a  -  radius  of  the  earth)1,  and  to 

consider  that  terrestrial  magnetic  field  in  the  first  approximation, 
it  carries  dipole  character,  so  that  at  the  height/altitude  $  above 


the  earth's  surface 
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where  tj  =  $/s;  //,<,0)  —  f ield  on  the  earth's  surface,  then  expression  for 
the  coefficients  a,  /3,  7,  6,  e  with  an  accuracy  to  the  values  of 
order  tj  they  will  take  the  form 


a  —  U{U*  —  Ya) — X(U* — /2Kn)-j-f)TjKa(fy — IX),  (6) 
p  =  2 XI  Ya  |  (C.OT  +  C,n)  -  (7C, rn  +  6C,n) ) ,  (7) 

f  =  -s2  [2fy2(a-.v)-2^  +  .v/^]-f2A‘f;(f;-,Y)+ 

+  AT*  |-l+(C1/»-*-C,«),|  - 

-  2rjc?[20/2(t/  — ^)-26,K^4-.Y/>ij  -  ^S'YlWU  -  XI*)- 
-2f{XYl\C1m{Cxm-^Cin)  +  Z\-\-[-{Cxm-\-C^\\%  ’  (8) 

8  =  -  2X1  Yl  |(5*  +  2 t,Ci)  (C,m  +  Ctn)  -  7,S2  (7C,m  -f  6C») ,  (9) 
e  =  (U  ~  X)  \U2 (54  +  4tAS0  -  2 XU  (52 -f  2t,C?)  +  AT*  j  + 

+  5s  {X Yl  1 1  —  (C,m  -f-  C»J  +  27iC,m  (C,m  +  C./i)  — 

-  6^(1  -  (Ctm  +  C2/i)'] )  -  UYl  (5 2  +  2i)C?  -  6t,S2)|  + 

4-2riC'fKa  (*(1  —  (CjW  -|-  C,/;)2)  -  iAS2),  (10) ' 

where 


52-1 -Cl;  CJ=C?  +  C|;  52  =  1  —  C2,  Yl  =  (^fji, 

~~  2 
w«a—  the  value  of  gyrof  requency  on  the  earth's  surface  x~~^r  ^ 

plasma  frequency;  ('=  1  +  / .  v9.m>— effective  collision  frequency;  l 

m,  n  -  direction  cosines  of  vector  with  the  axes  of  coordinates  p,  o 

z. 


FOOTNOTE  1 .  Let  us  note  that  in  work  [3]  the  coefficient  17,  which 
considers  the  final  curvature  of  the  ionosphere,  derived  from  the 
geometric  examination,  is  obtained  erroneously  (2$/a  instead  of  $/a) 
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ENDFOOTNOTE . 
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From  the  examination  of  equation  (4)  it  is  possible  to  obtain 
the  information  about  the  bias/displacement  of  the  roots  of  Booker's 
equation  [1],  caused  by  the  presence  of  curvature  in  the  ionosphere 
and  by  the  heterogeneity  of  terrestrial  magnetic  field.  Any  simple 
analytical  expressions  for  this  it  is  impossible  to  obtain,  but 
qualitatively  it  is  possible  to  say  that  the  values  of  the  roots  of 
equation  (4)  decrease  in  comparison  with  the  flat/plane  case,  but  the 
effect  of  the  final  curvature  of  ionosphere  and  magnetic  bump  of  the 
earth/ground  on  the  height/altitude  noticeably  is  manifested  in  the 
region,  where  approximation/approach  WKB  becomes  inapplicable. 

In  conclusion  I  express  deep  gratitude  to  the  docent  of  the 
department  of  radiophysics  of  LGU  V.  V.  Novikov  for  the  granting  of 
theme  and  direct  management/manual  with  the  execution  of  this  work. 
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